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ABSTRACT 

This  paper  presents  a  frequency-domain  technique  for  estimating 
distributed  lag  coefficients  (the  impulse-response  function)  when  observations 
are  randomly  missed.  The  technique  treats  stationary  processes  with  randomly 
missed  observations  as  amplitude-modulated  processes  and  estimates  the 
transfer  function  accordingly.  Estimates  of  the  lag  coefficients  are  obtaine 
by  taking  the  inverse  transform  of  the  estimated  transfer  function.  Results 
with  artifically  created  data  show  that  technique  performs  well  even  when  the 
probability  of  an  observation  being  missed  is  one-half  and  in  some  cases  when 
the  probability  is  as  low  a6  one-fifth.  The  approximate  asymptotic  variance 
of  the  estimator  is  also  calculated  in  the  paper. 
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Introduction: 

Least-squares  estimation  of  the  lag  coefficients  of  a  distributed  lag 
model  is  not  a  straight-forward  regression  problem  when  the  sample  has  missing 
observations.  Even  though  the  normal  equations  can  be  computed  using  sums  of 
the  available  cross  products,  the  estimators  in  the  solution  of  these 
equations  are  not  necessarily  unbiased.  In  any  event,  a  special  computer 
program  is  required  to  compute  the  normal  equations  of  a  distributed  lag 
linear  model  for  this  case. 

In  this  paper,  we  present  a  frequency  domain  approach  to  the  estimation 
of  lag  coefficients  when  observations  are  randomly  missed.  The  estimators  are 
easy  to  compute  if  a  discrete  Fourier  transform  program  is  available,  as  is 
the  case  for  most  time  series  software  libraries.  Estimates  of  the  lags 
computed  by  our  approach  have  common  large  sample  variance.  This  is  a  useful 
property  for  identifying  model  structure. 

Our  paper  proceeds  as  follows.  In  the  first  section  we  review  the 
^’-equency-domain  approach  to  the  estimation  of  distributed  lag  models.  In  the 
second,  we  discuss  the  estimation  of  spectra  and  distributed  lags  with 
missing  observations.  In  the  final  section  we  present  some  results  of  using 
our  suggested  technique  on  artifically  created  data. 


1 .  A  Frequency  Domain  Approach: 

Our  discussion  of  the  frequency  domain  approach  to  the  estimation  of 
distributed  lags  follows  the  commonly  used  non-parametrlc  approach  (Jenkins 
and  Watts,  1968,  Chapter  10).  Let  (x(t)}  and  (e(t)}  be  two  independent,  mean 
zero,  stationary  time  series,  and  let 

E[x(t+r)x(t)]  «  ox(T)  for  all  t 

and 

2 

°e  t  -  0 

E[e  (t+r )e(t)]  -  , 

0  T  *  0 

so  that  (e(t)}  Is  white  noise.  For  simplicity,  set  the  time  unit  equal  to  the 
sampling  interval,  so  that  t  takes  on  Integer  values.  The  time  series  (y(t)} 
is  related  to  the  other  two  series  according  to: 

00 

y(t)  -  l  h(k)x( t-k)  +  e(t),  (1) 

k  *-“> 

where  (h(k)}  Is  an  absolutely  summable  sequence  whose  transfer  function  Is 
zero  outside  the  band  -it  <  u>  <  it.  We  will  also  assume  that  h(0)  *  0t  so  that 
changes  In  x(t)  cannot  have  an  Immediate  impact  on  y(t). 

From  (1)  we  obtain 

Sxy(w)  “  H(u))Sx(o)),  -it  <  (D  <  it,  (2) 


where  Sx(w)  and  SXy(oj)  denote  the  own  spectrum  of  {x(t)>  and  the  cross- 


spectrum  between  (x(t)}  and  (y(t)},  respectively,  and 


H(u>)  =  \  h(k)  exp(-io)k)  (3) 

k=-oo 

is  the  transfer  function  of  (h(t)}.*  Upon  rearranging  (2), 


H(oj)  “  Sxy((o)/Sx(u). 


(4) 


Now,  suppose  that  we  observe  the  series  (y(t)}  and  {x(t)}  at  the  same 
times  t=0, 1, . . . ,n-l.  From  this  set  of  observations  we  can  obtain  estimates 
of  the  spectrum  of  {x(t)}  and  the  cross-spectrum  between  {x(t)}  and  {y(t)} 
as  follows.  Let  X(ujc)  denote  the  discrete  Fourier  transform  of  {x(t)}  for 
tie  observations  at  the  angular  frequency  *  2rrk/n;  i.e., 


X(u>k)  =>  \  x(t)exp(-ia)jct). 

t=*0 


(5) 


*Let  Ya(k)  be  the  absolutel  •  summable  covariance  function  of  the  mean 
zero,  stationary  time  series  {a(t)}.  We  use  the  convention  that 


Sa(<o)  “I  Ya(k)e_i“k, 

k— ” 


— IT<(jO<1T 


defines  the  spectrum  of  (a(t)}.  Under  this  convention 


Ya(k)  =  —  /  Sa(<u)e^u)kdaj. 
2  ir  _ 
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r 

a 


Thus,  we  can  estimate  Sx(wk)  as 
*  1  d 

Sx(<»k)  "  ~  1  IX((ujc+j)l2  »  (6) 

8  j— d 

where  s**2d+l.  Similarly,  if  we  let  Y( uik)  denote  the  discrete  Fourier 
transform  of  {y(t)}  at  the  angular  frequency  u>k,  we  can  estimate  as 

'  1  d 

SxyCw^)  ■  ~  \  X(wk+j)Y  (aik+j),  (7) 

s  j— d 

where  an  asterisk  denotes  complex  conjugate. 

Thus,  we  can  form  an  estimate  of  H(ui)  by  substituting  (6)  and  (7)  into 
(4)  to  obtain 

a  A  a 

H(u))  -  S^oO/S^aO.  (8) 


From  (10.3.14)  in  Jenkins  and  Watts  we  know  that  when  n  and  s  are  large 

*  1  -2 
E|H(o)k)  -  H(wk)l2  «  -  !H(u)k)|2[Yw(a>k)-l],  (9) 

s  ’ 

2 

where  Yxy^k^  lfl  tlle  squared  coherency  between  (y(t))  at  (x(t)}  at  frequency 
uk;  l.e., 

Y  (l*)k  )  *  I S  jjy  ( (iik  )  I  /Sjj  (ojj^  )Sy(cok) . 


An  estimate  of  (h(k)}  can  be  obtained  by  taking  the  discrete  inverse 
transform  of  H(u);  l.e.,  we  can  estimate  (h(k)}  as 


■»  1  n-1  * 

h(k)  -  -  1  H(w.)exp(iukj). 

n  j-0 

Given  the  smoothing  in  the  numerator  and  denominator  of  (8) 
-  1  m-1 

h(k)  »  -  l  H(u)48)exp(iu)4sk), 

™  j-0 


(10) 


(10') 


I 
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where  m  *  n/s.  If  n  +  «  as  n  +  »  ,  then  the  expected  value  and  variance  of 
the  right  hand  side  of  (10')  converge  to  the  expected  value  and  variance  of 

a 

(10).  Since  {H(tujs):  j=0,...,m-l}  are  asymptotically  Independent  as  m  ♦  «°, 
It  follows  from  (9)  that 

Var[h(k)]  «  —  l  E|H(ti) js)  -  H(tojs)|2  (11) 

m2  j-0 

-  ---  l  lH(u)j8)l2[Yxy  >j8)  "  1] 

nr  s  j*0 


- -  /  |H(u))l2[Yxy(ai  -  l]du. 

2itn  0 


The  estimators  of  (h(k)}  have  the  same  variance  at  all  lags  since  this 
variance  Is  independent  of  k.  Equation  (11)  suggests  that  the  variance 

A 

of  (h(k)}  can  be  estimated  as 


Var[h(k)]  -  -  \  |H(w .) | 2[yxv(w s )  “lj.  (12) 

n  j=*0  J 


m 
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2.  Missing  Observations: 

Our  method  of  estimating  (h(k)}  in  the  case  of  missing  observations 
proceeds  along  lines  similar  as  those  above.  However,  it  first  involves 
finding  a  method  of  estimating  Sx(ai)  <ind  ^xy  (a>)  for  the  case  of  randomly 
missing  observations.  Our  approach  is  to  follow  Bloomfield  (1970)  who  applies 
Parzen's  (1963)  approach  to  amplitude-modulated  stationary  processes  to  the 
estimation  of  single  series  spectra  in  the  case  of  randomly  missing 
observations. 

Consider  the  case  of  the  (x(t)}  series.  The  amplitude-modulated  series 
{x* (t)}  is  constructed  by  replacing  the  missing  observations  in  the  original 
series  by  zero,  the  mean  of  (x(t)}.  In  other  words,  we  define  the  amplitude 
modulating  series  {z(t)>  *  (x1 (t)/x(t)}  as 

1  if  x(t)  is  observed 

z(t)  -  (13) 

0  otherwise. 

Thus, 

x'(t)  *  z(t)x(t).  (14) 

Following  Bloomfield  and  Scheinok  (1965),  we  assume  that  the  pr >cess 
which  causes  observations  to  be  missing  is  a  stochastic  process  which  is 
independent  of  (x(t)}.  In  addition,  we  assume  that  (z(t)}  has  the  properties 
that  for  all  t, 

P[z(t)  »  1]  -  px 
p[z(t)  —  0 ]  —  1  —  Px 

and 


Cov[ z(t ) ,z(t  +  t)]  »  oz(t). 
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Thus  (z(t)}  is  weakly  stationary,  md 

E[z(t)z(t  +  t)1  “  <jz(t)  +  p2  for  all  t. 

1  x 

By  the  independence  of  {x(t)}  and  {z(t)}»  it  follows  that  E{x'(t)}  * 
E{x(t )>E{z(t)}  =0.  Consequently, 

E{x’(t)x'(t  +  t)}  -  E{x(t)x( t  +  t)}  E{z(t)z(t+  -)}  (15) 

-  ox(r)  [p2  +  oz(t)] 

X 

It  then  follows  that  the  spectrum  of  {x'(t)J  is 

Sx»(m)  ■  p2Sx(w)  +  (2x)-^  J  Sz  (oj— to  *  )  Sx(u')diD'.  (16) 

X  “IT 

Similarly,  we  assume  that  there  is  an  amplitude-modulating  series  ?(t) 
such  that 

1  if  y(t)  is  observed 

C(t)  -  <17> 

0  otherwise 

which  gives  rise  to  the  amplitude-modulated  series 

y ' (t)  -  ? (t)y(t).  (18^ 


We  assume  that  (?(t)}  is  independent  of  both  (x(t)}  and  {y(t)>  and  that  the 
series  has  the  properties  that  for  all  t, 


P[c(t)  -  l]  -  Py 
P[c(t)  -  0]  -  1  -  Py 


Cov[c(t).z(t  +  t)]  -  oz5(t). 


Thus, 


E[c(t)z(t  +  t)]  -  02?(t)  +  PxPy* 

Since  (x(t)}  and  {y(t)}  are  independent  of  both  (z(t)}  and  (?(t)},  it  follows 
that 


E{x'(t)y'(t  +  t)}  -  E{x(t)y(t  +  t)}  E{z(t)C(t  +  t)}  (19) 

"  Cfxy(t)[pxPy  +  OZC(T)]. 

It  then  follows  that  the  cross  spectrum  of  (x' (t)}  and  (y’(t)}  is 

i  ,‘n 

Sx*yi((o)  *  PxPysxy^u^  +  (2w)-1  /  S2^  (uj-o)' )SXy(a)' )du' .  (20) 

Thus,  using  estimators  similar  to  (6)  and  (7)  we  could  estimate  Sxi(w), 

Sx*y»((o),  S2(w),  and  (oi)  from  the  available  data  and  use  (16)  and  (20)  to 
solve  for  Sx(w)  and  S-gyiu).  These  estimates  could  then  be  substituted  into 
(8)  and  (h(k)}  estimated  by  taking  the  inverse  transform  of  the  result,  as  in 
the  case  when  no  observations  were  missing.  Note  that  our  analysis  does  not 
require  either  that  (x(t)}  and  (y(t)}  be  observed  at  the  same  times  or  that 
the  same  random  process  generates  the  missing  observations  of  the  two  series. 

In  addition,  there  is  no  requirement  that  the  random  processes  governing 
the  observability  of  the  underlying  series  be  independent,  so  that  the 
analysis  does  not  preclude  the  case  in  which  (x(t)}  and  (y(t)}  mist  be 
observed  at  the  same  times. 

In  general,  (16)  and  (20)  can  be  solved  for  Sx(u>)  and  Sj^w)  using  the 
method  given  by  Bloomfield.  However,  in  the  present  paper  we  will  only 
discuss  the  special  case  in  which  the  series  (z(t)}  and  (c(t)}  are  white 
noise.  In  this  case,  (16)  and  (20)  can  be  considerably  simplified  and  the 
solutions  for  Sx(o>)  and  SXy(oi)  become  much  easier. 

V 

.  M 


When  the  binomial  distributions  generating  the  missing  observations  are 


white  noise. 


Px(1-Px)  T  =■  0 

0  (t  )  * 

z  0  otherwise 


and 


az c  T  “  0 

ozc(t)  - 

0  otherwise 


(21) 


(22) 


Under  these  assumptions  (16)  and  (20)  become 


and 


Sx»  (a) ) 


p2sx(a>)  +  Px(l“Px)0^ 
x  x 


(23) 


Sx*y  i  (w)  m  PxPysxy  (2A) 

respectively.  Thus,  we  find  that  In  the  case  in  which  (z(t)}  and  (c(t)}  are 
white  noise,  the  spectra  of  the  amplitude-modulated  series  are  linear 
transformations  of  the  original  series.  Substituting  (23)  and  (24)  into  (4) 
yields 

2 

H(oi)  ”  (Px^Py )  {[sx»y»(d))  “Pz^^xy  ] /[ Sx i  (lo )  _  Px^^Px^x) (23) 
Thus,  if  we  estimate  Sxi(u)k)  as 

Sx.(wk)  -  -  1  lx*(uk+1)!2  (6’) 

8  j— d 


in 


and  Sxty  i  (iok)  as 

Sx,y«(iok)  -  -  l  X' (a)k+j)V'*(a)k+j),  (7') 

s  j— d 

we  can  estimate  H((uk)  as 

a.  a  a  a  a  a  a  A  a  A  ^ 

H(uk)  *  (px/py){  [Sx'y'^k^  “  azi;axy  ] /[  Sx' ^k^Px^  1_Px>°x]  > •  <26) 

Further,  if  oz)-  *  0  or  -  0,  It  Is  shown  in  the  Appendix  that 

^  2 

E  |H(a)k)  -  H(iok)  | 2  *  -  |H(u)k)  !  2{yx'y  •  (oik)  +  [(  “]^  "!}•  (27) 

s  px  SxUk) 

This  approximation  is  good  when  p/s  >  >  1,  where  p  *  min(px,Py). 

Once  again  an  estimate  of  {h(k)}  can  be  obtained  by  taking  the  discrete 

a 

inverse  transform  of  H(w)  as  given  by  (26).  Further,  we  could  obtain  the 

A 

approximate  variance  of  (h(k)}  by  repeating  the  steps  used  to  obtain  (11) 
using  (27)  in  place  of  (9).  This  variance  is 

2  1  ° 

Var[h(k)]  =  -  /  |H (u) i 2{yx« y * (m)  +  [( - )  ~  “]2  ”  l}dw.  (28) 

2irn  o  Px  sx'Mk' 

It  can  be  estimated  in  the  same  way  that  (11)  was  estimated  above. 
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Once  again,  we  note  that  the  asymptotic  variance  of  h(k)  Is  the  same 
for  all  lag  coefficients.  Further,  note  that  if  ax^/S(oj)  and  ay2/Sx(<o)  are 
close  to  unity  for  all  w  and  px  =  py  =  p,  the  variance  In  the  case  missing 
observations  will  be  more  than  times  that  for  the  case  In  which  no 
observations  are  missing. *  Further,  when  p  gets  small,  the  standard  errors  of 
{h(k)}  are  of  the  order  (p/n)~*. 


*lf  cr2£  =  0  or  Oxy  =  0,  it  is  easy  to  show  tha- 

2  2  r  (1_Px>  Ox2  ,r  d-Py)  °y2  , 

Yx'y'(“)  "  + - “K1  + - T7"^}  *  (29) 

SX(W)  Py  Sy\&J/ 


2  2 

Thus,  if  Sx(u))  -  ax  and  Sy(ui)  «  Oy,  the  terras  in  brackets  in  (29)  are 

-1  -1 

approximately  equal  to  px  and  Py  ,  respectively.  Thus,  for  this  for  this 

2  2 
case  Yx'y'Cm)  =  (pxPy)Yxy 


3.  Results  with  Artlf ically  Created  Data: 

In  order  to  evaluate  our  frequency-domain  approach  to  estimating 
distributed  lags  when  observations  are  randomly  missing,  we  used  It  on  some 
artifically  created  data.  Specifically,  we  generated  {y(t)}  according  to  (1) 

using  the  weights 

0.25  k  -  1,  5 

0.75  k  -  2,  4 

h(k)  *  (30) 

1.0  k  =  3 

0.0  otherwise. 

The  errors  (e(t)}  were  computed  by  a  normal  N(0,1)  pseudo-random  number 
generator.  For  all  experiments  the  number  of  observations  was  12000,  which 
Is  approximately  the  number  of  days  since  1947.  The  Independent  Bernoulli 
trials  used  to  create  (z(t)}  and  (c(t)}  were  obtained  using  a  uniform  (0,1) 
pseudo-random  number  generator.  In  the  case  of  (z(t)},  if  the  value  of  the 
t-th  random  number  exceeded  px,  then  z(t)  was  set  equal  to  zero.  Otherwise, 
z(t)  =  1.  A  similar  procedure  was  used  to  generate  (c(t)}  with  py  used  in 
place  of  px.  In  all  experiments  (z(t)}  and  {^(t)}  were  independent.^  All 
spectral  estimates  were  obtained  by  smoothing  pertodograms.  The  weights  used 

A  A 

in  the  smoothing  of  S x * y « ( (d^ )  and  Sxi(u)k)  are  discussed  below.  In  the 

^The  specific  normal  pseudo-random  number  generator  used  was  the  FTGEN 
subroutine  in  the  International  Mathematics  and  Statistics  Library  (IMSL). 

The  specific  uniform  (0,1)  pseudo-random  number  generator  used  was  the  GGUBS 
subroutine  of  IMSL. 


calculation  of  (26)  actual  value?;  of  px,  py,  and  oz^  were  used  instead  of 
estimated  ones.^  The  variance  of  {x(t)}  was  estimated  using  only  the  non-zero 
values  of  {x' (t)}. 

The  first  set  of  experiments  were  performed  letting  (x(t)}  be  white  noise 
normal  random  variables  with  unit  variance.  The  (x(t)}  were  generated  with 
the  same  computer  subroutine  used  to  generate  (e(t)}.  The  results  for  each  of 
up  to  ten  runs  are  presented  in  Table  1.  The  values  used  for  px  and  py  are 
given  in  the  column  labeled  "Percentage  Observation. ”  In  all  runs  a  999  point 
moving  average  was  used  to  smooth  the  periodogram  to  obtain  the  spectral 
estimates.  The  row  labeled  "Squared  Error  Fit"  contains  the  values  of 
11999  - 

l  [h(k)  -  h(k)]2, 
k-0 

which  is  the  sum  of  squared  differences  between  the  actual  and  estimated  lag 
coefficients. 

There  are  three  major  points  concerning  these  results.  The  first  is  that 
the  estimated  lag  coefficients  track  the  actual  lag  coefficients  quite  well  as 
long  as  px  and  py  are  0.2  or  above.  This  can  be  seen  both  in  the  low  values 
of  the  squared  error  of  fit  and  the  high  degree  of  correspondence  between  the 
means  of  the  estimated  lag  coefficients  and  their  actual  values. 

The  second  point  is  that  the  variances  of  the  parameter  estimates 
increase  as  the  percentage  of  missing  observations  increasesas  (28)  and  (29) 
indicates  they  should.  This  can  be  seen  in  the  increasing  standard  errors  of 
all  of  the  estimates  as  the  "Percentage  Observation"  decreases.  Further,  it 

^Since  the  values  which  would  have  been  obtained  if  px,py  and  ozC  had 
been  estimated  were  close  to  the  true  values,  the  results  would  not  change 
much  if  estimated  values  of  those  parameters  had  been  used  Instead. 
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is  interesting  to  note  as  long  as  px  and  py  are  0.2  or  above,  these  standard 
errors  are  roughly  constant  cross  coefficients  as  (27)  indicates  they  should 
be. 

The  third  point  is  that  there  is  a  downward  bias  in  the  estimates  of  h(l) 
to  h(5),  the  non-zero  lag  coefficients.  Our  conjecture  is  that  this  downward 
bias  is  caused  by  the  bias  imparted  to  the  spectral  estimates  due  to  including 
999  terms  in  the  moving  average  used  to  smooth  the  periodogram.  To  test  this 
conjecture,  we  redid  several  of  the  runs  with  Sxtyi(i»>^c)  and  Sx(<j>k)  estimated 
using  only  99  terms  in  the  moving  average.  These  results  are  contained  in 
Table  2.  As  can  be  seen,  this  modification  does  decrease  the  downward  bias 
for  the  cases  px  =*  Py  ’  .9  and  px  =*  Py  *  .5.  However,  it  is  also  apparent 
that  for  this  case  the  estimated  lag  coefficients  track  the  true  lag 
distributions  well  only  when  no  more  than  50  percent  of  the  observations  are 
missing.  Thus,  we  have  the  classic  trade-off  between  bias  and  variance  due  to 
the  fact  that  increasing  the  number  of  terms  used  to  smooth  the  periodogram 
decreases  the  variance  but  Increases  the  bias  of  the  estimates  of  the 
spectrum. 

We  also  evaluated  our  procedure  when  (x(t)}  was  colored.  Specifically, 
using  the  same  nNormal  N(0,1)  pseudo-random  number  generator  as  before  we 
generated 


x(t)  -  px(t-l)  +  v ( t ) , 


(19) 


2 

with  p  «  0,5  and  crv  ■  1.  We  also  generated  enough  observations  on  x(t)  prior 


to  those  used  to  estimate  the  lag  coefficients,  so  that  the  choice  of  initial 
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x(t)  should  not  Influence  the  results.  Our  results  for  this  case  are 
presented  In  Table  3.  Once  again  we  find  that  we  are  able  to  track  the  lag 
d  strlbutlon  well  as  long  as  no  more  than  50  percent  of  the  observations  are 
missing.  Further,  the  results  for  this  case  continue  to  show  that  variance  of 
the  estimates  increase  as  px  and  py  decrease  and  that  including  more  points  in 
the  moving  average  used  to  smooth  the  periodograms  increases  the  bias,  but 


decreases  the  variance  of  the  estimated  coefficients. 
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APPENDIX 


The  Approximation  (27)  for  the  Mean  Squared  Error  of  H (to ) 

A  A 

Let  ex*y»(<d)  »  Sxty«(u))  -  Sx«yi(w)  and  ex«(a))  ■  Sxi(u>)  -  Sxt(io).  Assume 
that  the  cumulants  of  (x(t))  and  {y(t)}  satisfy  condition  (4.3.10)  in 
BriLlinger  (1975). ^  From  Theorems  4.3.2  and  4.4.1, 


Eex'y'Go)  *  0,  Eex«(u))  «  0, 

E|cxiyi(u))l2  *  -  Sxi(a))Syi(ai), 

Eex»2(w)  *  -  Sx«2(w), 

8 


(Al) 

(A2) 


(A3) 


and 

*  1  * 

Eex»y«  (u)6x*  (u)  5  *  Sx*yi(u))Sx*(u))  ( A4 ) 

3  s  3 

for  large  n  (the  approximations  are  of  order  n“^).  Since  the  variances  of 

A  A  A 

ex«y»  and  ex»  are  of  order  s”^,  whereas  the  variances  of  px>  py,  and  ox2  are 
of  order  n”*,  we  can  substitute  px,  py,  and  ox2  for  their  estimates  in  (26) 
without  affecting  the  order  of  the  approximation  when  n  »  s. 

Since  0£2  ■  0  or  Oxy  *  0, 

A 

H(w)  »  (px/py)[Sx«y»  (u>)  +  €x»y>(u)]/[Sx»((o)  +  e  x»  (w)  -  c]  (A5) 


^This  condition  holds  if  (x(t)}  ind  (c(t)}  are  Guassian  ARMA  processes 
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where  c  =*  px(l  -  Px)ax^*  The  first  order  term  In  a  Taylor  series 
approximation  of  (A5)  Is 

A 

H(w)  -  H(u))  «  pxpv_1(Sx*  (w)  -  c)"1  [ex'v*(tt») 

(A6) 

-  (Sxi(a))  -  c)-1Sxiyt(i»))ex«  (co > ] 

It  follows  from  (23),  (24),  (A2),  and  (A3)  that  this  approximation  is  good  If 
1/p/s  is  small,  where  p  =  min(px,  Py).  Applying  (A2)  and  (A3)  to  (A6), 

A 

E|H(w)  ~  H  (oj  )  I  ^  ~  s-1  [pXPy“^  (Sx<  (u>)  “  c)"2] 

[sx«  (w)Sy  »(<ij)  “2  |Sx*y  *  (a;)  I  2(SX>  (u)  -  c)-1Sxt(u>) 

2 

+  |Sxiy  i  (u)  |  2(SX»  (id)  -  c)-2Sxi(w)j 

-2 

3  s”1[PxPy~l(sx,(<*))  “  c)"2lsx>y'(u)l  ]2  [yx.y.(u) 

-2(SX»(<»))  -  c)"^Sx*(u))  +  (Sx»  -  c)~2  Sx<(<i))] 

-2 

=  s-1|H(<o)|2{Yx'y»(i*>)  -  1  +  [i  -  (Sx«(ai)  -  c)_ISxt  (aj)]2 
-2 

=  s~MH(a))|2{Yxty»(aj)  -  1  +  [c(Sx(to)  -  c)-1]2} 

-2  1-PX  ax2 

■  8-1  |H(u)|  2{Yx*y*  (a>)  -  1  +  [( - ) - ]2}. 

Px  Sx(at) 
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